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Abstract 



In the present paper we address the interaction of pairs of charmed mesons with hidden 
charm in a finite box. We use the interaction from a recent model based on heavy quark 
spin symmetry that predicts molecules of hidden charm in the infinite volume. The energy 
levels in the box are generated within this model, and from them some synthetic data are 
C^ • generated. These data are then employed to study the inverse problem of getting the 

energies of the bound states and phase shifts for DD or D*D* . Different strategies are 
investigated using the lowest two levels for different values of the box size, carrying a 
^ ■ study of the errors produced. Starting from the upper level, fits to the synthetic data are 

carried out to determine the scattering length and effective range plus the binding energy 
of the ground state. A similar strategy using the effective range formula is considered with 
^ ' a simultaneous fit to the two levels, one above and the other one below threshold. This 

^f^ ■ method turns out to be more efficient than the other one. Finally, a method based on 

tJ" , the fit to the data by means of a potential and a loop function conveniently regularized, 

turns out to be very efficient and allows to produce accurate results in the infinite volume 
^r . starting from levels of the box with errors far larger than the uncertainties obtained in 

jyv , the final results. A regularization method based on Gaussian wave functions turns out 

to be rather efficient in the analysis and as a byproduct a practical and fast method to 
calculate the Liischer function with high precision is presented. 

X 

1 Introduction 

The determination of the hadron spectrum from lattice QCD (LQCD) calculations is attracting 
many efforts and one can get an overview on the different methods used and results in the recent 
review [1] . One of the tools becoming gradually more used is the analysis of lattice levels in terms 
of the Liischer method [2,3]. This method converts binding energies of a hadron-hadron system 
in the finite box into phase shifts of hadron-hadron interaction from levels above threshold, or 
binding energies from levels below threshold [4-6]. From the phase shifts one can get resonance 
properties, and there are several works that have recently applied these techniques to study the 
p resonance [7-15]. There exist other resonances far more difficult to get with this approach like 
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the ai(1260), which was also attempted in [14] (see also its determination using finite volume 
calculations with effective field theory in [16]). Scalar mesons have also been searched for with 
this method [17-20] and gradually some calculations are being performed for systems in the 
charm sector [21-27]. From another field theoretical perspective, finite volume calculations have 
also been devoted to this sector in [27-29]. In [28] the KD, rjDs interaction is studied in finite 
volume with the aim of learning about the nature of the D*o(2317) resonance from lattice data.^ 
The theoretical model used in [28] is taken from [31], where the 15*0(2317) resonance appears 
dynamically generated from the interaction of KD, rjDs and other less relevant channels. In this 
latter work, a scalar hidden charm state coming from the DD interaction with other coupled 
channels was also found, which qualifies basically as a DD quasi-bound state (decaying into 
pairs of lighter pseudoscalars) . Although not reported experimentally, support for this state 
has been found in [32] from the analysis of the data of the e"'"e~ — )■ J/ipDD reaction of [33]. 
From the effective field theory point of view, this state has also been reported in [34,35], using 
light SU(3)-fiavour and heavy quark spin symmetries to describe charmed meson-antimeson 
interactions. 

The purpose of the present paper is to study the interaction of DD and D*D* using a field 
theoretical approach in finite volume in order to evaluate energy levels in the finite box which 
might be compared with future LQCD calculations. The paper also presents a strategy to better 
analyze future lattice results in order to get the best information possible about bound states 
and phase shifts in the infinite volume case from these lattice data. For this purpose we shall 
use the model of [35], although most of the results and the basic conclusions are independent 
of which model is used. 

As to the method to obtain the finite volume levels and the inverse problem of obtaining 
the results in the real world, phase shifts and binding energies, we shall follow the method 
of [36] where a reformulation of Liischer approach is done based on the on shell factorization 
of the scattering matrix that one uses in the chiral unitary approach [37-42]. This method is 
conceptually and technically very easy and introduces improvements for the case of relativistic 
particles (although we shall not make use of the relativistic version in the present paper). Some 
works using this formalism can be found in Refs. [16,28,29,43-47]. 



2 Formalism: infinite volume 

In this section, we briefiy review the formalism of Refs. [34,35], where an effective field the- 
ory incorporating SU(3)-light fiavour symmetry and heavy quark spin symmetry (HQSS) is 
formulated, to study charmed meson-antimeson (generically denoted here HH' , with H^ H' = 
D,D*,Ds, D*) bound states. The lowest order (LO) contribution of the interaction is given 
by contact terms, and the symmetries reduce the number of independent low energy constants 
(LECs) of the approach to only four. Other effects, like one-pion exchange or coupled chan- 
nel dynamics, are shown to be sub-leading corrections to this order. Still, coupled channels 
will be considered explicitly when the mass difference between the thresholds is not negligi- 
ble compared with the binding energy of the molecules considered. To fix the four constants 
of the approach, one assumes the molecular nature of some XYZ states, namely, X(3872), 
X(3915) and y(4140). The fourth input of the model is the isospin violating branching ratio 



^The first hint, though indirect, of the nature of the D*q(2317) as a mostly KD bound state from lattice 
data was presented in [30]. There, lattice calculations of the scalar form factors in semileptonic pseudoscalar- 
to-pseudoscalar decays were used to extract information about the corresponding elastic S'-wave scattering 
channels. 



of the decays X(3872) — > J /'^u and X(3872) — > J j'^p (for a different approach to this issue 
and the X(3872) — )■ J/'^^ decay see [48]). For further details on the formahsm we refer to 
Refs. [34,35,49]. We will adapt here the formalism to a more adequate (for the problem in 
hands) T-matrix language. 

Since we are dealing with heavy mesons, we use a non-relativistic formalism. In our normal- 
ization, the S'-matrix^ for an elastic HH' scattering process reads 

S{E) = e2^^(^) = 1 - t^T{E) , (2.1) 

where the modulus of the momentum k = \k\ is given by k"^ = 2p{E — mi — m2), and fi is the 
reduced mass of the system of two particles with masses mi and m2. In Eq. (2.1), 6 is the phase 
shift, and we can write: 

T = -^sin5e^'^ , (2.2) 

fik 

T-i = -^cot5 + z^. (2.3) 

The expression for the T-matrix is given by: 

T"\E) = V~\E) - G{E) , (2.4) 

with V the potential (two-particle irreducible amplitude) and G a one-loop two-point func- 
tion. This equation stems from a once-subtracted dispersive representation of T~^{E) (see for 
instance Sec. 6 of Ref. [39]), or equivalently, from the N/D method [50] equations, when the 
left-hand cut is neglected or included perturbatively [40,51-53]. The loop function G provides 
the right-hand cut and the contribution of the left-hand cut should be included in the potential 
V. As mentioned above, we will follow here the approach of Refs. [34,35], and we will approxi- 
mate V by its LO contribution in the l/rriQ expansion (with rriQ the mass of the heavy quark). 
Thus, we are completely neglecting the left-hand cut. 

The loop function G needs to be regularized in some way. Typical approaches are once- 
subtracted dispersion relations and sharp cutoffs. Here, instead, we are following the approach 
of Refs. [34,35], in which the loop function is regularized with a Gaussian regulator. For an 
arbitrary energy E, we find 

f dV g-2(q2-fc2)/A2 

^ ^ ~ J (27r)3 E- mi -7712- q^f^p + iO+ 

-^e-M%i|.(y2./A)-.g. (2.5) 

with (j){x) given by: 

0(x) = Te^'dy . (2.6) 

Jo 

Note that, the wave number /c is a multivalued function of E, with a branch point at threshold 

[E = m,i + m,2)- The principal argument of [E — mi— 7712) should be taken in the range [0, 2ti[. 

The function k(l){\/2k/A) does not present any discontinuity for real E above threshold, and 



^We will always consider S*— wave meson-antimeson interactions, and thus the spin of the molecule will always 
coincide with the total spin of the meson-antimeson pair. The partial waves ^■^+^Lj are then '^^'^^Sj^s- For 
simplicity in what follows, we will drop all references to the L, 5" and J quantum numbers, both in the S and 
T matrices. 



G{E) becomes a multivalued function because of the ik term. Indeed, G{E) has two Riemann 
sheets. In the first one, ^ Arg(£' — ?72i — 7712) < 27r, we find a discontinuity G i{E + ie) — G i{E — 
ie) = 2ilTaGj{E + ie) for E > {mi + 7712)- It guaranties that the T— matrix fulfills the optical 



theorem. For real values of E and below threshold, we have k = i J—2jj,{E — rrii — 7712). Poles 
below threshold in the first sheet correspond to bound states. In the second Riemann sheet, 
27r ^ ATg{E — mi — 7712) < 47r, we trivially find Gii{E — ie) = Gi{E + ie), for real energies and 
above threshold. 

The Gaussian form factor enters Eq. (2.5) in a way that is unity for on-shell momenta, and 
hence the optical theorem ImT~^ = /ifc/(27r) is automatically fulfilled. A Gaussian regulator is 
also used in Ref. [54] to study the A resonance in finite volume. 

The cutoff A is a parameter of the approach, and, hence, the theory depends on it. This 
dependence, however, is partially reabsorbed in the counter-terms of the theory, as long as 
one chooses a reasonable value for it, not beyond the high-energy scale of the effective field 
theory [55-58]. 

In the approach of Refs. [34,35], the potential V is taken as 

V{E) = e-2'='/^'c(A) , (2.7) 

where C is the proper combination of the four different counter-terms for each considered chan- 
nel HH'. Explicit expressions can be found in Appendix A. The dependence of the counter- 
term on the ultraviolet (UV) cutoff A should cancel that of the loop function G, such that 
G{Eb)V{Eb) becomes independent of A, when Eb is the energy of the bound state used to de- 
termine the counter-term.^ For other energies, there will exist a remaining, unwanted/unphysical, 
dependence of the T matrix on the cutoff. This is due to the truncation of the perturbative 
expansion (see discussion in Ref. [34]). Up to this point, we have discussed only the case of 
uncoupled channels, but the generalization to coupled channels is straightforward.'^ 
Finally, above threshold the effective range expansion reads: 

kcot6 = — + -rA;2 + • • • , (2.8) 

where a and r are, respectively, the scattering length and the effective range. From Eqs. (2.3) 
and (2.4) we can calculate the theoretical predictions for these effective range parameters, 
obtaining: 

87r / 1 ziA \ , , 



3 Formalism: finite volume 

In this section, we follow the steps of Ref. [36] to write the amplitude in a finite box of size 
L with periodic boundary conditions, denoted by T. Since the potential does not depend on 



■^The mass of a bound state is thus given by T ^{Eb) = for Eb < rrii + m2, which is the equivalent of 
Eq. (10) in Ref. [35]. 

'*One just has to rewrite the T- matrix as T = (I — VG)^^V , where V and G are now matrices in the coupled 
channels space. 



L, one only has to replace the loop function G with its finite volume version, G, in which the 
integral over momentum q is replaced by a discrete sum over the allowed momenta, 

f-^{E) = V-\E) - G{E) , (3.1) 

1 -2{g2-fc2)/A2 

G(E) = -^ y ^rr- , (3.2) 

^ ' L^y E -nil -7712- g2/2/i 

where the (quantized) momentum is given by: 

27r 
^= fi^ neZ^. (3.3) 

Ij 

Now, the energy levels in the box are given by the poles of the T-matrix, V~^ = G. For the 
energies of these levels in the box, the amplitude in the infinite volume is recovered as: 

T~\E) = V~\E) - G{E) = G{E) - G{E) = 5G{E) . (3.4) 

Since the G function is regularized (either in the box or in the infinite volume) with a Gaussian 
regulator, the difference above depends explicitly on the cutoff A. This remaining non-physical 
dependence on A quickly disappears as the volume increases. Indeed, we find that it is exponen- 
tially suppressed and that it dies off as exp (— L^A^/8) (see Appendix B). Thus, it is clear that 
in this context, we can end up the renormalization program just by sending the UV cutoff to 
infinity. This will allow to obtain the physical T— matrix, independent of any renormalization 
scale, for the energy levels found in the lattice Monte Carlo simulation (finite box). 

For the practical calculations that we will show in what follows, the A dependence is already 
negligible when A > 1 GeV even for the smallest volumes considered in this work (the limit 
A — ;■ oo is effectively achieved for such values). In this limit, Eq. (3.4) becomes the Liischer 
equation [2,3], as we discuss in certain detail in Appendix B. The results of the appendix 
also show that the inclusion of a Gaussian regulator is a quite efficient technique, from the 
computational point of view, to evaluate the Liischer function Zoo(l,fc^) used in [3]. Finally, 
from Eqs. (2.3) and (3.4), we can write: 

kcoi5 = -— lim Re (G{E) - G{E)) . (3.5) 

y(i A-i-OO V / 



4 Results 

We present in this section the results obtained with the formalism outlined in the previous 
section. We first discuss the results obtained by putting the model of Refs. [34, 35] directly in 
the box. That is, we study the volume dependence of the molecules found in Ref. [35], predicting 
thus the existence of sub-threshold levels (asymptotically different of threshold) for the different 
channels, which have a clear correspondence with the hidden charm molecules reported in [35]. 
This is done in Subsec. 4.1. 

Our purpose in Subsecs. 4.2, 4.3 and 4.4 is to simulate a realistic situation in a LQCD study, 
where one would obtain different energy levels (one or two) for different sizes, L, of the box. 
To do so, we generate "synthetic data" from the exact levels that we obtain from the model 
of Refs. [34,35]. We take five different values of Li, in the range L77i.,^ = 1.5 to 3.5. From the 
calculated levels, we obtain randomly shifted levels (in a range of 5 MeV), and assign an error 
of 10 MeV to each of these points, except the last one for which we have assumed an error of 



only 8 MeV to prevent it from crossing the threshold. Next we use a Monte Carlo simulation, 
to estimate the errors on the determination of observables (the phase shifts, for instance) when 
the energy levels are obtained with a certain statistical error. Specifically, we study in these 
subsections the 1 = J^'" = 0"*""^ DD channel. 

In Subsec. 4.2, the Liischer formalism to study the phase shifts calculated from Eq. (3.5) is 
applied to the synthetic levels above threshold that we find for the different studied channels. 
From these phase shifts, we calculate the effective range expansion parameters, and use them to 
determine the masses of the bound states. In Subsec. 4.3 we adopt another strategy to extract 
information from the generated levels. Namely, we consider a potential which parameters are 
then fitted to reproduce the synthetic levels (above and below threshold, simultaneously). With 
this potential, we can make predictions in the infinite volume case, and thus we end up with 
another determination of the masses of the predicted bound states. We shall see that this 
method allows one to obtain better results (better central value and smaller errors) for the mass 
of the bound state than the previous one. We then analyze in detail which are the differences 
of both approaches. In Subsec. 4.4 another method is proposed, in which the effective range 
approximation is retained for the inverse of the T-matrix amplitude, but fitting directly the 
energy levels instead of the phase shifts, and studying simultaneously the levels above and 
below threshold. In this case, then, we notice that the precision achieved for the mass of the 
bound state is similar to that obtained with the potential analysis. 

In Subsec. 4.5, we analyze in a more quantitative way, the qualitative arguments given in 
Subsec. 4.1, where the behavior of the sub-threshold levels is discussed. We offer a method to 
discriminate between those levels that produce bound states in the limit L -^ oo and those 
that do not, and hence tend to threshold in the infinite volume limit. This method allows the 
extraction of the mass and the coupling of the bound state in the infinite volume limit. 

All these methods are applied in Subsec. 4.6 to the bound state present in the / = 
J^'-' = 2++ D*D* channel. The difference with respect to the case used as an example in the 
previous subsections is that the state is now weakly bound (the binding energy is only around 
2-3 MeV), so that we can compare how the methods exposed above work for this case. 

4.1 The model of Ref. [35] in a finite box 

In Figs. 1 and 2 we present the dependence of the energy levels on the size of the finite box, 
as calculated from Eqs. (3.1) and (3.2), for the different channels studied in Ref. [35]. We 
have fixed the potential in the different channels by means of the central values given in this 
reference for the various counter-terms, and collected here in Appendix A, Eqs. (A.3)-(A.6). 
When needed, we have also implemented in the finite box a coupled channel formalism. The 
solid lines correspond to the case A = 1 GeV, whereas the dot-dashed lines to A = 0.5 GeV. For 
comparison, we also show, with the horizontal dashed lines, the involved threshold energies. 
We just show those energy levels that can be identified with bound states {k"^ < 0) in the 
infinite volume case. That is, their asymptotic L — ^ oo value approaches the bound energies 
given in [35], and thus they are different from threshold. Of course, one has this latter piece 
of information from the calculations of the model in an infinite volume, but this would not be 
the case in a lattice simulation. Let us focus, for simplicity, in the 1 = case, shown in the left 
panel of Fig. 1. The large L asymptotic behavior can be well appreciated in some cases like the 
O'^"*" DD or 1"* D*D molecules. However in other cases, it might be difficult to discriminate 
between a real bound state and a threshold level, even for quite large values of the box size L. 
Clear examples are the 1^"*" D*D or the 2+"'' D*D* molecules (similar examples can be found in 
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Figure 1: Volume dependence of the / = (left) and / = 1/2 (right) molecules predicted in 
Ref. [35]. The horizontal dashed lines show the different thresholds involved (when the charge is 
not explicitly given, we are displaying the thresholds associated to the different charge channels). 
The solid lines correspond to the levels found in the box for A = 1 GeV, whereas the dot-dashed 
ones stand for those obtained with A = 0.5 GeV. Over the right axis we mark with arrows the 
masses of the bound states as predicted in the infinite volume case and A = 1 GeV. The J^*^ 
quantum numbers of the different channels are indicated beside the arrows. 



the different isospin-strangeness channels), which in the infinite volume case are loosely bound. ^ 
Thus, we see a well known result from Quantum Mechanics; the smaller the binding energies, 
the larger become the L values needed to reach the asymptotic behavior. From this study, we 
conclude that in a lattice simulation when dealing with states that are at least bound by some 
tens of MeV, one might safely discriminate them by using box sizes of the order of Lm.^- ~ 3. 
However, the lattice study of other states less bound (binding energies below 4 or 5 MeV, as in 
the case of the X(3872) or the 2+"*" resonances) might require significantly larger volumes. To 
achieve more accurate predictions for the former and solve the problem for the latter ones, we 
follow different approaches in the following subsections. 

Finally, we note that some of the levels in Figs. 1 and 2 are not realistic, in the sense that 
they would mix with other levels generated by channels with the same quantum numbers, but 
lower thresholds. That is the case, for example in / = 0, of the 0++ D*D* at -E ~ 3920 MeV, 
that would mix with some higher levels of the DD channel. Indeed, it is to be expected that 
these bound states would acquire some width due to the coupled channel dynamics. Still, it is 
possible that these states could appear as more or less stable energy levels.^ 



^The first one corresponds to the X(3872) resonance that has been observed close to the D^D^* threshold [59] 
(see also a recent determination and discussion of other experiments in Ref. [60]) and it has been a hot topic 
for both the experimental and theory communities since its discovery. The 2+"'' state is a HQSS partner of the 
X(3872) molecule which dynamics, at LO in the heavy quark expansion, is being determined by precisely the 
same combination of counter-terms that appear in the X(3872) channel. Given the discovery of the AT (3872) 
resonance, the existence of the 2++ state, either as a bound state or a resonance, is therefore a quite robust 
consequence of HQSS [34,35]. 

^In Ref. [61] the D*D*, DID* states arc studied with the interaction taken from the extrapolation of the 
local hidden gauge approach to the charm sector, which also respects HQSS. The coupling to DD and DgDg is 
allowed and generates a width of about 50 MeV for the most bound state, the one with / = and J^'-^ = 2++. 
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Figure 2: The same as in Fig. 1 for the 1=1 (left) and hidden strangeness (right) molecules 
predicted in Ref. [35]. 



4.2 Inverse analysis: phase shifts 
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Figure 3: Some energy levels for the / = 0, J^'^ = 0++ DD interaction as a function of the box size 
L. The levels obtained with the model of Refs. [34,35] in a box for A = 1 GeV are shown with (red) 
solid lines, while the generated levels for some particular values of L (synthetic data points, see 
the text for details), together with their assigned errors are displayed with black circles. The non- 
interacting energies (mi + 771,2 + (27r/L)^n^/2/i with 7i = 0, 1) are shown with (blue) dash-dotted 
lines. The error bands around the solid lines are obtained from the fit to a potential discussed in 
Subsec. 4.3. They have been obtained by considering pairs of fitted parameters (1/Coa, A) that 
provide values of x^ that differ from the minimum one by less than one unit (x^ ^ Xmin + !)■ 



We start discussing the case of the isoscalar 0++ DD interaction. Some levels found from 
the model of Ref. [35] in a finite box/ obtained as the zeros of Eq. (3.1), are shown with a 
(red) solid line in Fig. 3. The synthetic levels generated from them and our choice for their 
errors are shown with points. Recall that we give an error of ±10 MeV to these points trying 
to simulate a realistic situation in a LQCD study, where these levels will be determined with 



^In what follows we will use an UV cutoff A = 1 GeV when presenting results deduced from the model of 
Ref. [35], both for finite boxes and in the infinite volume case. Other cutoffs compatible with the effective theory 
designed in [34, 35] give rise to similar results. 
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Figure 4: Phase shifts obtained for / = J^*^ = 0++ DD interaction. The points stand for 
the phase shifts calculated from the synthetic energy levels displayed in Fig. 3 using Eq. (3.5). 
The green dashed line and its associated error band corresponds to the effective range analysis of 
Subsec. 4.2, while the blue solid line and its error band stand for the results obtained by fitting 
a potential discussed in Subsec. 4.3. The phase shifts in the infinite volume are very similar to 
the latter ones, so we do not show them. In both cases, the error bands have been obtained by 
considering pairs of fitted parameters [(l/a,r) for the effective range fit and (1/Coa, A) for the 
case of the potential fit] that provide values of x^ that differ from the minimum one by less than 
one unit (points included in the dark blue y^ ^ Xmin + 1 ellipse displayed in Fig. 5). 



some statistical uncertainties. From the upper level and the Liischer's formula, Eq. (3.5), we 
find the phase shifts shown with points in Fig. 4. The errors in the phase shifts in this figure are 
determined by recalculating them, through Eq. (3.5), with different values of the upper level 
energy, E, randomly taken within the error intervals displayed for each of the synthetic data 
points in Fig. 3. 

We could also obtain the scattering length and the effective range parameters either from 
the determined phase shifts, or from Eqs. (2.8) and (3.5). Actually, combining these two latter 
equations we have. 



ReSGi = lim Re (G(E) - G(E)) = -— ( -- + -rk^ + 
A-i>oo V / /i V a 2 



(4.1) 



for the upper energy levels, E, determined in finite boxes of different sizes. We have obtained 
1/a and r from a x^— linear fit to the five data points generated for lie 6Gl using the five 
synthetic upper energy levels^ shown in Fig. 3. We find 



1 



0.62 ±0.25 im-\ 



0.53 ±0.18 fm 



(4.2) 



with a linear Gaussian correlation coefficient R = 0.83. From the above result, we find 



(4.3) 



These values are to be compared with those obtained in the infinite volume model, Eqs. (2.9) 
and (2.10), with parameter G = Coa(A = 1 GeV) = —1.024 fm^, which turn out to be: 



Oth = 1-38 fm , Tth = 0.52 fm . 



(4.4) 



®To estimate the errors in Re SGl for each of the synthetic energy levels considered, we follow a procedure 
similar to that outlined above for the phase shifts. Thus, we let the synthetic energy level vary within the error 
interval displayed in Fig. 3 and find the range of variation of Re SGl ■ 
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Our fitted values are compatible with the theoretical ones, but have sizeable errors, although 
the correlation is large. Performing a standard analytical continuation of Eqs. (2.3) and (2.8) 
below the DD threshold, we estimate the position of the X(3715) bound state, 

E = 3721+^° MeV , (4.5) 

whereas the value found in Ref. [35] is E = 3715li5 MeV.^ The binding energy, 5 < 0, is 
obtained from Eqs. (2.3) and (2.8), upon changing k — )• in, and imposing T~^ = 0, 



k2 



B=-, K= V '_ (4_g) 

2/i r 

To estimate the uncertainties in Eq. (4.5), we have performed a Monte Carlo simulation taking 
into account the existing statistical correlations between 1/a and r. We quote a 68% confident 
interval (CL), but with some caveats as we explain next. Note that 2r/a is not far from unity 
and within errors it can be even bigger, which means that we can get some events in the Monte 
Carlo runs (around 25%) with 1 — 2r/a < 0, for which we set the square root to zero. Thus, 
the lower error quoted in Eq. (4.5) is somehow uncertain, since the above procedure tends to 
accumulate events around 3695 MeV. On the other hand, for the cases with 1 — 2r/a > 0, but 
small, the two roots of k in Eq. (4.6) are not so different, and hence there is some ambiguity 
in the binding energy B (we choose the smallest value of k). Note that, although the value of 
E obtained with its errors seems quite accurate, when one considers it relative to the binding 
energy B, we find a large dispersion, since the DD threshold is at around 3734 MeV. 

Finally, if we decrease the error of the synthetic energy levels from 10 MeV to 5 MeV, then 
the errors of the phase shifts as well as those of the threshold parameters are also reduced 
approximately to half of their previous values, and the predicted mass is more accurate, E = 
3723!tii MeV (and now only for around 6% of the Monte Carlo events, 1— 2r/a become negative). 
This should give an idea of the precision needed in the determination of the energy levels in 
order to have an appropriate determination of the mass. 

In the next sections we discuss different alternatives that allow to achieve a better precision. 

4.3 Inverse analysis: fit to a potential 

We now consider another approach to analyze/use the synthetic levels that we generated in the 
previous section. Here again we aim to falsifying real data obtained from LQCD Monte Carlo 
simulations for various finite volumes. The analysis of phase shifts in the previous subsection 
necessarily takes into account only the level above threshold in Fig. 3. It is then convenient 
to develop an approach that could simultaneously make use of all available levels. Thus, we 
propose to describe all levels using a potential, Eqs. (2.7) and (3.1), fitting its parameters for 
such purpose. We adopt here an approach where we fit a counter-term C = Coa defining the 
potential^'' and the UV cutoff A (involved in the finite box loop function and in the potential 
(see Eq. (2.7)) to the synthetic energy levels shown in Fig. 3. Thus, the x^ function is then 
given by: 



^The errors calculated for finite volume quantities in tliis work refer to the statistical uncertainties we generate 
in the synthetic data. The errors quoted from Ref. [35] refer instead to the uncertainties in the determination 
of the constants appearing in the potential. 

^"We follow here the notation of [34, 35] where the counter-term that appears in this channel is called Coa 
(see Appendix A). 
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where E^^ir {Li 



are the first two energy levels calculated from the HQSS potential, with pa- 



^thV 

rameters C^a and A, in a finite box of size Lj. On the other hand, E'^'"''' and ^El"'''' are the 
synthetic levels, that we have generated, together with their assigned errors. Here, the super- 
script j = 0, 1 refers to the two levels shown in Fig. 3. The fit parameters, 1/Coa and A, obtained 
in the best fit are 



1 



a 



-0.93 ± 0.20 fm-^ A = 970 ± 130 MeV 



Oa 



(4.8) 



with a linear Gaussian correlation coefficient R = —0.98. These errors, and the correlation 
coefficient, are calculated from the hessian of x^ at the minimum. However, since the fit is not 
linear, these errors are slightly different from those obtained requiring x^ ^ Xmm + -'-• This latter 
requirement gives the following non-symmetrical errors: 



1 

Coa 



-0.93l°i? fm-^ 



A = 970t}^[! MeV 



(4.9) 



From Eq. (4.8), we find C^a = — 1.08J;o;29 fm . The central values of both, the counter-term and 
the UV cutoff, agree well with those of the original model of Ref. [35], Coa = —1.024 fm and 
A = 1 GeV, used to generate the synthetic levels. However, as expected, the two parameters are 
strongly correlated. This is further discussed in Appendix C A contour plot of the x^ function 
in the (1/Coa, A)— plane is shown in Fig. 5, that manifestly shows the correlation. 

On the other hand, the fitted parameters of Eq. (4.8) predict a value for the mass of the 
bound state oi E = 3715^6 MeV (68% CL, obtained from a Monte Carlo Gaussian simula- 
tion keeping the statistical correlations) in the infinite volume case. The central value agrees 
remarkably well with the value obtained from the model of Ref. [35], E = 3715 MeV, and cer- 
tainly much better than that obtained with the phase shift analysis carried out in the previous 
subsection {E = 3721l';;25 MeV). The errors found now are also significantly smaller. 

Finally, we have calculated error bands for the predicted finite box levels and phase shifts by 
quantifying the variations that are produced in these observables when one randomly considers 
pairs (1/Coa, A) of parameters that provide x^ ^ Xmin + 1 (points included in the dark blue 
ellipse displayed in Fig. 5). These error bands are shown in Figs. 3 and 4, respectively. 




600 



-1.6 



-1.4 



-1.2 



-0., 



-0.6 



-0.4 



Figure 5: Contour plots for the x^ function defined in Eq. (4.7) for the / = J^'^ = 0++ DD 
channel. The dashed line shows the correlation predicted from the model of Ref. [35], Eq. (C.2). 



The circle represents the central value taken for that model, Coa 
while the square stands for the results of the best fit of Eq. (4.8). 



1.024 fm^ and A = 1 GeV, 
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4.4 Inverse analysis: effective range 

We have seen in Subsecs. 4.2 and 4.3 that the fit of the synthetic energy levels with a potential 
leads to better results for the mass of the bound state than those obtained from the fit to the 
phase shifts (deduced from the upper level) with the effective range expansion. We believe that 
the reasons for this improvement are basically three. First, the potential fit takes into account 
both levels, above and below threshold, while the phase shifts analysis takes into account just 
the upper level. Second, in the potential fit, the "observables" are the energy levels, while, in the 
phase shifts analysis, the quantity that enters in the x^ function is k cot 6, and the propagation of 
errors can lead thus to worse determinations of the parameters. Third, the analytical structure 
of the inverse of the amplitude is different in both approaches. In the effective range approach, 
ones truncates a series up to k"^, while in the potential fit one is effectively including further 
terms beyond the latter ones. Indeed, the full loop function G is taken into account. To study 
the importance of the first two points, we follow here another approach, in which we shall keep 
the effective range approximation for the amplitude, but fit the energy levels (above and below 
threshold) instead of the phase shifts obtained from the above threshold level. 

The effective range expansion for the inverse of the T-matrix is written from Eqs. (2.3) and 
(2.8): 



f -7^ ( -- + ^/c^ - ik] , k = \ 2ii(E - mi - mo) , E > mi + mo 
I -2^ [-a - i7 + 7j , 7 = V2^('Tii + m2-E) , E < mi + m2 

Now, the energy levels in the box are found for given values of a and r, by means of Eq. (3.4). 
It is to say, by numerically solving T~^{E) = SGl{E), similarly as it is done in the case of the 
potential, but now using Eq. (4.10) to model the T— matrix both above and below threshold. 
We will denote the levels obtained in this manner as £^thEF- To determine the best values of a 
and r, we consider thus a x^ function as Eq. (4.7), where the E{'' are still the synthetic levels 
we have generated, but replacing the E^^y by -EthEF) calculated as explained above. The values 
of the best parameters are: 

- = 0.70 ± 0.07 fm-\ r = 0.56 ± 0.07 fm , (4.11) 

a 

with a linear Gaussian correlation coefficient R = —0.6.^^ Hence, we obtain: 



a = lASroil fm . (4.12) 

The errors calculated in this way are clearly smaller than those displayed in Eq. (4.2) with the 
phase shifts analysis carried out in Subsec. 4.2. Also, the central value of the scattering length 
agrees better with the theoretical one, Eq. (4.4). These improvements have a clear impact in 
the determination of the mass of the bound state, which is now Eb = 3716]'^5 MeV (68% CL), 
with smaller errors and better central value than those obtained with the phase shifts analysis 
in Eq. (4.5). In Fig. 6 we show a comparison of the ellipses in the (l/a,r) parameter space 
determined by the condition x^ ^ Xmin + 1 ^^ ^^e fits of Eqs. (4.2) and (4.11). There, it can be 
clearly appreciated the significant improvement achieved by fitting directly to both, lower and 
upper energy levels, instead of fitting to the phase shifts deduced from the latter levels. Finally, 
we must point out that the determination of the energy levels obtained with this method are 
very similar to those obtained in Subsec. 4.3 by introducing of a potential, and shown in Fig. 3. 



In this case, the errors calculated from the requirement x ^ Xmin + 1 almost coincided with those given 
here in Eq. (4.11). 
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Figure 6: Comparison of the determination of the effective range parameters, 1/a and r, with 
the methods explained in Subsecs. 4.2 (red dashed hne) and 4.4 (solid blue line). The ellipses 
are obtained from the condition y^ ^ Xmin + 1 in each case. The central values of each fit are 
represented with points. Wc also show, for comparison, the theoretical values of the parameters, 
as given in Eq. (4.4), with a black square. 

Actually, the differences between the upper and lower energy level curves (and their error bands) 
deduced from both methods would not be easily appreciated in Fig. 3. For this reason, we have 
not shown in this figure the results obtained from the method discussed in this subsection. 



4.5 Inverse analysis: bound state fit 

We have discussed in Subsec. 4.1 the volume dependence of the sub-threshold levels that arise 
when we put the model in a finite box. For cases with V^ < 0, the potential is attractive, and 
hence, a bound state in the infinite volume case may arise. Whether it is bound or not in the 
infinite volume case, there would appear a sub-threshold level for finite volumes. It was argued 
in Subsec. 4.1 that it may be not very clear, at first sight, if these levels tend to the threshold 
energy or to a bound state in the L — )■ cx) limit. To circumvent this problem, we suggest here 
a method to study this volume dependence. By subtracting Eqs. (2.4) and (3.1), we can write 
the amplitude in the finite box as [62]: 



T-i = T-^ - 6Gl, 6Gl = lim 6G = lim {G-G 

A— 5>oo A— s-oo 



(4.13) 



A bound state with mass Eb appears as a pole in the T-matrix, thus in the vicinity of the pole, 
we can approximate: 



TiE ~ Eb) 



9' 



(4.14) 



E-Eb 

where the ellipsis denote regular terms in the Laurent series of the amplitude. The coupling 
can also be calculated analytically. 



g^ = lim (E - Eb)T{E) or 



1 dT-i(^) 



E^Eb 



dE 



(4.15) 



E=Eb 



The volume dependence of the sub-threshold level in the finite box, given by the equation 
T~^{E) = is then dominated by this bound state, and hence: 



f~\E) 



E-E 



B 



5Gl{E) = 



(4.16) 
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from where one can write: 

E{L) = Eb + g^5GL \E{L), L] . (4.17) 

This equation is a reformulation of a similar result obtained in Refs. [4-6].^^ The coupling g 
obtained here is related to Z^ of Ref. [6]. Note, however, that Eq. (4.17) is appropriate as long 
as Eq. (4.14) is sufficiently accurate to describe the infinite volume T— matrix for the energy 
levels found in the lattice simulation (i.e., energies for which T~^ vanishes). Hence, the larger 
the box sizes, the better Eq. (4.17) will perform.^^ We extract the mass and the coupling of the 
bound state from a fit to the sub-threshold level in Fig. 3 with the following x^ function. 



5 ITT^fTA _ 77^(0)^2 



E(Li)-E , 
where E{L) is given by Eq. (4.17).^^ The best x^ is Xmin — O-^? and the parameters so obtained 



are: 



Eb = 3712 ± 6 MeV, g'^ = {2.8 ±2.1) GeV-\ (4.19) 



to be compared with those obtained with the model in the infinite volume case, Eb = 3715 MeV 
and g^ = 2.6 GeV"\ 

It could be that, for the case of weakly bound states, the error bars on the energies overlap 
with the threshold and it is difficult to determine if one has a bound state or not. A weak 
attractive potential that does not bind in the infinite volume case, still provides a level below 
threshold for finite volumes, and the energies go to threshold as L -^ oo. For this case we 
proceed as follows. The volume dependence of this level would be given by: 

f^ f '^ , ^12 , \ 2/i 1 2 






L^k"^ 2na \ An 

with some coefficients a and P, disregarding exponentially suppressed terms. In the above 
equation, we have explicitly separated the threshold singularity of 6Gl, 

1 1 2/i 



L^ E — nil — ""^2 L'^k'^ 



^^In the same line as in those references, but using boosted reference frames, in Ref. [63] hnear combinations 
of energy levels are suggested to reduce the volume dependence. Our method does not rely upon the analytical 
form of the volume dependence to make extrapolations since for every L considered, the exact L dependence is 
provided by JGl (£(£), L). 

^^On the other hand, for very small binding energies, some subtleties appear, because the coupling g^ tends 
to zero as the mass of the bound state approaches the threshold [53,64,65]. We will discuss this issue at length 
in Subsec. 4.6. 

-'^^It is worth noting the following technical detail. In principle, E{L) should be extracted for each Li as the 
implicit solution in Eq. (4.17) for given Eb and g^. For practical purposes, though, it is more convenient to 
obtain E(L) by plugging into the right-hand side of this equation the values of E^ and Li that we arc fitting 
to. If Eq. (4.14) is accurate enough, both methods are equivalent, as long as the effects in Eq. (4.17) of the 
statistical fluctuations of the measured lattice levels are sufficiently small. In that case, the results for Eb and 
g^ should not be very different, as we have checked. Indeed, the best fit results given in Eq. (4.19) have been 
obtained within this approximation. However, this approximation cannot be safely used when the bound state 
is placed close to threshold, because then 5Gl rapidly changes and statistical fluctuations in the determined 
lattice energy levels induce large variations in the right-hand side of Eq. (4.17). 
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Hence, the general behavior of this level would be: 

,2 2/i 1 



L^ A + Bk"^ 



(4.20) 



Now, this expression could be used in a x^ function as in Eq. (4.18) (with k^ < Q for the level 
below threshold). We have checked that, if we try to fit the energies of the lower level of Fig. 3 
with this formula, we get a much worse x^ value, discarding the possibility that there is not a 
bound state in the infinite volume. 

We have seen then that the method outlined in this subsection allows for a safe discrimi- 
nation between those levels that correspond to bound states in the infinite volume and those 
that do not, and it also gives a precise determination of the mass. It is also worth noting that 
the errors for the mass of the bound state are similar to those obtained with the analysis from 
a potential in Subsec. 4.3, and smaller than those calculated with the phase shift analysis in 
Subsec. 4.2. 

4.6 The case of the / = J^^ = 2++ channel 

We now repeat the same analyses carried out above but now for the case of the bound state 
present in the 1 = J^^ = 2~^~^ D*D* channel. As already mentioned, the 2"''+ state is a HQSS 
partner of the X(3872) molecule which dynamics, at LO in the heavy quark expansion, is being 
determined by precisely the same combination of counter-terms that appear in the X(3872) 
channel. The existence of the 2~^~^ state, either as a bound state or a resonance, is a quite robust 
consequence of HQSS [34,35], and it can be certainly subject to experimental detection. It is 
also worth to discuss this channel in this context because, contrary to the case analyzed before, 
we have here a very weakly bound state. In the calculation of Ref. [35], and also in the results 
of Subsec. 4.1, shown in the left panel of Fig. 1, charged and neutral coupled channels were 
studied, because the gap between both thresholds is indeed larger than the binding energy of 
the bound state. The mass of the bound state is Eb = 4013.2 MeV, whereas the neutral and 
charged thresholds are located at 4014.0 MeV and 4020.6 MeV, respectively. Here, however, 
in order to discuss the problem in simpler terms, we will consider only an uncoupled channel 
problem with / = 0, and use isospin average masses, keeping the relevant counter-term, Cq 
in the nomenclature of Ref. [35], to the same value, namely, Cq = —0.731 fm^. In this way, 
the threshold is located at 4017.3 MeV, whereas the mass of the bound state now becomes 
Eb = 4014.6 MeV. The first two energy levels obtained with this simplified model are shown 
with (red) solid lines in Fig. 7. As before, for the following statistical analyses we consider 
the synthetic levels shown in this figure with points. The centroid of these points is randomly 
shifted in the range ±5 MeV, and they are given an error of ±10 MeV. With the points of the 
upper level, we generate the phase shifts shown in Fig. 8 with points, through Eq. (3.5). We 
can now obtain the scattering length and effective range as in Subsec. 4.2, that turn out to be: 

- = 0.41 ± 0.30 fm-\ r = 0.67 ±0.19 fm , (4.21) 

(Ji 

with a linear Gaussian correlation R = 0.81. From the above fitted value for 1/a, we find: 

a = 2.4^2.4 fj^ ^gg^^ (^L) , (4.22) 

while the theoretical values, obtained from Eqs.(2.9) and (2.10), are: 

ath = 3.0 fm , rth = 0.58 fm , (4.23) 
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Figure 7: Same as in Fig. 3, but for the / = 2++ D*D* interaction. 
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Figure 8: Phase shifts for the / = J^'^ ~ 2++ D*D* channel interaction. The points stand for 
the synthetic phase shifts generated from the upper energy level of Fig. 7. The (green) dashed line 
correspond to the effective range fit, while the (blue) solid line correspond to the potential fit. The 
associated error bands are obtained by considering randomly chosen pair of parameters [(l/a,r) 
and (1/Co, A) for the effective range and potential fits, respectively] that satisfy x^ ^ Xmin + 1- 



which agree with the above determinations, within errors. The phase shifts obtained with these 
parameters, and the associated error bands, are shown in Fig. 8 with (green) dashed hnes, and 
they satisfactorily reproduce the synthetic data. The mass of the bound state turns out to be 
then Eb = 4013_"i8 MeV. Recall that the caveats raised in Subsec. 4.2 apply here. 

The next step is to consider a potential fit, as in Subsec. 4.3, where the free parameters are 
Co and A. From a best fit, as described in Subsec. 4.3, we find: 



1 
C'o 



-1.4 ±0.5 fm" 



A = 1020 ± 240 MeV , 



(4.24) 



with a linear Gaussian correlation coefficient R = —0.97. The non-symmetrical errors given by 
the condition x^ ^ Xmin + 1 ^^^ instead: 



Cn 



-1 -4-0:7 fm- 



A 



10201?^!^ MeV . 



(4.25) 



The energy levels obtained with these parameters are shown in Fig. 7 with a (black) dashed line, 
although they are so similar to those of the exact model that they mostly overlap. Also shown 
in the figure are the error bands generated from the errors of the parameters, as described 
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previously in Fig. 3. With this potential, we can also calculate the phase shifts, which are 
shown in Fig. 8 with (blue) solid lines, and also the associated error band. The quality of both 
descriptions, that of the effective range and that of the potential, are very similar, and indeed 
both lines are very similar to the one of the exact model, and hence we do not show the latter. 
The value of Co deduced from Eq. (4.24) is Cq = -O.Tltoig fm^ (68% CL), in good agreement 
with the one of the infinite volume model, Cq = —0.73. Finally, the mass of the bound state is 
given by Eb = 4014.3t5;4 MeV. We must stress again that the errors obtained with this method 
are smaller than those obtained with the phase shifts analysis. 

Now, we consider the method of Subsec. 4.4, in which the effective range expansion is used 
to study the levels below and above threshold. In this case, the best fit values that we obtain 
for the scattering length and effective range parameters are: 

- = 0.35 ±0.15 fm-\ r = 0.64 ±0.15 fm , (4.26) 

a 

with a linear Gaussian correlation coefficient R = 0.3. The non-symmetric errors stemming 
from the condition x^ ^ Xmin + 1 ^^^i^ o^^ to be: 

1 = 0.35+°i? fm-i , r = 0.64+°;^^ fm . (4.27) 

Propagating the correlated Gaussian errors of Eq. (4.26), we find: 

a = 2.9+g:^ fm . (4.28) 

We note that here, as it also occurred for the / = 0, J^^ = 2"'""'" case, the central values obtained 
with this method agree better with the theoretical ones, Eq. (4.23), than those obtained with 
the phase shift description, Eq. (4.21), and have smaller errors than the latter ones. The mass 
of the bound state obtained is Eb = 4014.2l4;g MeV, which is better determined than that 
obtained by means of the phase shifts analysis, and very similar to the one obtained with the 
potential method. 

Finally, we should proceed now with the analysis performed in Subsec. 4.5. However, there 
is a major difference in this case, namely, that the bound state is very close to the threshold 
in this case. It is known that, in this case, the coupling of the state tends to zero [53,64,65], 
and so additional terms in the Laurent series, Eq. (4.14), are relevant for energies not very far 
from the bound state mass.^^ Further, since we are considering an error of ±10 MeV in the 
energy levels, and we are trying to reproduce a bound state of binding energy of 2 — 3 MeV, we 
should expect a greater qualitative impact on the lower energy levels, which are the only ones 
considered in this method (note that the other methods examined above use always the upper 
levels as well). These considerations explain why, when performing such analysis, we obtain 
very bad results for the bound state mass and the coupling. Hence, we must conclude that this 
method can only be applied safely in the case of bound states that are not very loosely bound. 
On the other hand, including a background term would increase the number of free parameters, 
and so, the errors stemming from a best fit would be even larger. At this stage, the approaches 
in Subsecs. 4.2 to 4.4 would be more useful for cases like this one, in which the bound state is 
very close to threshold. 



-'^^As an example, consider a background term in the amplitude in Eq. (4.14), so that T = g'^ /{E—Eb)+I3+- ■ ■ . 
From the theoretical model, one can calculate g^ = 0.58 • 10^'^ MeV^^ and /3 — —0.68 • 10^^ MeV~^. For energies 
E ~ 3990 MeV, as we find for the lower level in Fig. 7 for Lrrijr = 1.5, we have \g'^/{E - Eb)\ < |/3|. 
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5 Conclusions 

In this paper we have addressed the interaction of heavy charmed mesons in the hidden charm 
sector where several bound states are produced using an interaction that is based on heavy 
quark spin symmetry. The interaction is then studied in a finite box and the levels expected 
from a lattice QCD calculation are evaluated for the DD, D*D* states and their SU(3) partners. 
Then the inverse problem is faced, generating "synthetic data" from the levels obtained and 
using different procedures to obtain the relevant magnitudes in the infinite space, phase shifts 
and binding energies for the bound states. Particular emphasis is done in the error analysis to 
establish the accuracy of the different methods. We use two levels for different values of the 
box size L, one below threshold and the closest one above threshold. One strategy is to use the 
Liischer formula to get phase shifts for each energy of the level above threshold. Another strategy 
is to use the effective range approximation, but fitting the scattering length and effective range 
to both levels (above and below threshold). The two methods work, but the latter one gives 
better determinations of the parameters (scattering length and effective range), but also of the 
mass of the bound state. Yet, the method that proves most efficient ^^ is to parameterize a 
potential and a regularizing UV cutoff for the meson meson loops and carry a fit to the data. 
Once the potential and the UV cutoff are determined one can evaluate the phase shifts and 
binding energies with much better precision than the one assumed in the "synthetic data". The 
UV cutoff is not needed if one considers levels of only one energy, but it appears when different 
energies are used in the fit, yet, we show that it is highly correlated with the potential. We also 
devoted particular attention to the case of weakly bound states, where special care must be 
taken. Finally, as a byproduct we present an efficient method to obtain the Liischer function, 
supported by an analytical study that allows one to truncate the sum by means of a Gaussian 
form factor and estimate the error induced by the truncation. 
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A Potentials for the HH interaction 

In this Appendix, for completeness, we briefiy review the formalism of the effective field theory 
derived in Refs. [34, 35] to study charmed meson-antimeson bound states. This effective field 
theory incorporates SU(3)-light fiavour symmetry and HQSS. In this context and at LO, the 
potential is given by contact terms, related to four independent LECs, namely Coa, Cob, Cia 
and Cihi in the notation of Refs. [34,35]. Thus, the potentials in the different channels will be 
given by different linear combinations of these LECs. The fit of these LECs to four experimental 



^^As one moves far away from threshold, any method based on the effective range approximation becomes 
less appropriated. The potential fit method improves on the effective range expansion since it includes the full 
loop function and not only the imaginary part fixed by unitarity. Thus, this latter method should work better 
when energies significantly lower or higher than the threshold are considered. 
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Table 1: Potentials for the HH' interaction for the 7 = (left) and / = 1/2 (right) different J^'^ 
channels. 
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Table 2: Potentials for the HH' interaction for the 7=1 (left) and hidden strangeness sector 
(right) different J^*^ channels. 

data, as stated in Sec. 2, allows one to fix the numerical value of the different counterterms. In 
the following, we summarize the form of the potentials for the different / and J^*" channels. 

In the isoscalar channel, the only involved LECs are Coa and Cob, and the way they appear 
in the different J^^ channels is shown in Table 1 (left panel). Potentials in the / = 1/2 and 
7 = 1 sector are the same, except in the channel where coupled channels must be considered 
as discussed below. These potentials only depend on Cia and Cib. The explicit expressions for 
/ = 1/2 are given in Table 1 (right panel), whereas those for 1=1 are shown in Table 2 (left 
panel). Finally, in the hidden strangeness sector, the four LECs appear, and the final potential 
is the arithmetic mean of the corresponding isoscalar and isovectorial interaction, as can be 
seen in Table 2 (right panel). Note that, for this table, we have defined Coia = \ {Coa + Cia) 
and Coib = ^ (C'o6 + Cib). 

However, coupled channels must be taken into account in two cases. The first one is that in 
which the mass difference of the charged and neutral channels thresholds is not negligible as 
compared to the binding energy of the state. The second case occurs when charge conjugation 
is not a good quantum number. The first scenario is important in the study of the DD* system 
with J^*^ = l^"*" and D*D* with J^^ = 2+"*". In this case, the potential will be the interaction 
between the charged and neutral channels and will be given by the 2x2 matrix: 



Vo 



IfCo + Ci Co-Ci 
2 \ Co — Ci Co + Ci 



(A.l) 



being Co = Coa + Cob and Ci = Ci^ + Cib. The second scenario where coupled channels are 
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important is because of the mixing between the DgD* and DD^ channels. In this case, the 
potential is: 

^- - ( % t: ) ■ f^'^' 

Using these combinations of counter-terms to describe the four input data, as explained in 
Sec. 2, the numerical values of the LECs for the two values of the cutoff considered in this work 
are the following: 

Co. = -3.366^°:°?^ fm^ (-1-024+^03 fm^) , (A.3) 

Co. = +1.673;°;°J2 f^2 (+0.293+°:°°^ fm^) , (A.4) 

C,^ = -1.76+°:i fm^ (-0.684+n^^ fm^) , (A.5) 

Ci, = +l.68t'oil fm' (+0.311+°:°! fm^) , (A.6) 

for A = 0.5 GeV (1 GeV). These are the values used through this work. The uncertainties in the 
above equations account for possible HQSS violations and errors in the input used to fix the 
counter-terms (see a detailed discussion in Ref. [35]). For simplicity in this exploratory work, 
we have ignored them. 



B Gaussian regulator and relation to Liischer formula 

In this Appendix, we discuss the details of Eq. (3.5) within a Gaussian regularization scheme. 
We also study the dependence of the function SG{E), that appeared in Eq. (3.4), on the UV 
cutoff A. For convenience, we re-write SG{E) as 



6G{E;A) = G{E)-G{E) 



SGa 






d'q 


g-2(q2-fc2)/A2 _ I 


{2tiYE- 

dV 


- nil — ^^2 — 2 — '" ^0"*" 

1 


{2^YE- 


- mi - m2 - 1^ + iO+ 



L3 ^^ E-mi-m2-i^ 
+ I^^i?-mi-m2-g"/^^-'^- ^ "^ ^^~^ ^^-^^ 

SGl 

The function 6G explicitly depends on the cutoff A, and this dependence is carried by the 
SGa term. On the other hand, the term 6Gl is well defined, and it is related to the Liischer 
function [36] (see discussion below). In the strict A — )■ oo limit, only the second term survives, 
which justifies our approach in Sec. 3. Still, for practical purposes, the limit A — )■ oo can only 
be achieved by taking A large enough, and then it is necessary to study the dependence of 
6G with A. Let us note that SGa has no poles, and hence it is exponentially suppressed with 
L according to the regular summation theorem [2,3]. For k^ > 0, E > nii + ni2, SGl is not 
exponentially suppressed for L — > oo and, in this case, SGl clearly dominates over SGa- 

However for k"^ < 0, SGl is also exponentially suppressed as L increases, and therefore one 
needs to explicitly calculate the dependence of SGa on AL. 

Let us calculate the derivative of SG with respect to A. Only SGa depends on A, and this 
latter function does it through the exponential function exp[— 2(g^ — /c^)/A^]. The derivative 
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brings down a factor (g^ — k'^) that cancels out the denominators. This greatly simplifies the 
calculation of both the sum and the integral. The latter one is trivial and it only amounts to 
the integration of a Gaussian function, while the former one, up to constant factors, now reads 



L3 



E' 



-2g2/A^ 



J2 e-2(^) xJ 



^3(0, e"^ 



(B.2) 



where we have used q'^ = ql + Qy + (fz ^^^ that the exponential of a sum is the product of 
the exponentials. This latter property allows to relate the sum in three dimensions to the cube 
of the sum in one dimension. In Eq. (B.2), 9^{u,a) is a Jacobi elliptic theta function [66]. It 
satisfies [67], 

^3(0,6--') 1 



^3(0,6-^/^') X ■ 



This allows us to write then: 

d 5G 



dA 



(27r)3/2 



[^ 2feVA2 



^3(0, e 



-A2l2/8\ 



(B.3) 



(B.4) 



We note that this equation is exact. The above equation converges rapidly to zero as the 
Gaussian cutoff increases, which shows that the limit A — )• 00 is effectively quickly achieved. 
To proceed further, we note that: 



[^3(0, a)]^ = l + 6a + 12a^ + 



00 

m=0 



and the coefficients Cm are nothing but the multiplicities of m = n^, n G Z^. Since a = e 
we can find the leading term in Eq. (B.4), 

dSG 



dA 



6fi (2k^ A^L^\ , ,^, 



(B.5) 



-a^lVs 



(B.6) 



Given that 5G = SGl for A — ^ cxd, we find keeping just the leading term: 

A'^Ln 



SG{E;A) =6Gl{E) + 



SGl{E) + 



Qj2 



(27r)3/2 






(B.7) 



3/i 
2^ 



e erfc 



' AL 

2V2 



+ i 



V2k' 



-ikL 



erfc 



f AL .^/2k' 



and, then, its asymptotic behavior is: 



5G{E;A) =5Gl{E) 



24/i e 



(27r)3/2 AL2 



2 r2 



2(A;2L 



L2A2 



O (A"^) 



fB.8l 



where 0{A~^) refers to (k/A)^, [k"^ / L"^) / A'^ and 1/(LA)^, and the ellipsis stands for terms that 
are more exponentially suppressed (the next one would take the form e~^ ^ Z^). Given the form 
of the L suppression, the Gaussian regularization scheme does not introduce any spurious terms 
that would dominate over^^ the physical contribution 6Gl, as long as A is sufficiently large. 
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As already mentioned, for fc^ > SGl is not exponentially suppressed for L 
expect SGl to decrease as exp(— |fc|L). 



while for fc^ < 0, we 
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Indeed, one can efficiently compute the regularized Liischer function by means of the Gaussian 
regulated 5G loop function. The Liischer function^^ is related to the loop functions by means 
of [36]: 

Thus, for a mildly large value of A, SGl{E) can be approximated by the Gaussian regulated 
6G{E, A) function, up to corrections suppressed by the exponential factor e s~ (see Eq. (B.8)), 

4^Zoo(l, P) = -^^ (sGiE^A) + O (e-^)) (B.IO) 

which in turn provides -Zoo(l, k'^) with enough accuracy in a computationally easy way. 



C Cutoff effects and relation to dispersion relations 

In this Appendix and to better frame the approach followed in Subsec. 4.3, the existing correla- 
tion between the constant of the potential and the cutoff is addressed in detail. We also discuss 
the relation of our approach to other approaches in which the loop function is calculated from 
a dispersion relation. 

We recall Eqs. (2.5) and (2.7) to expand the inverse of the amplitude in powers of k"^. For 
more general purposes, we consider a potential in which the factor that multiplies the Gaussian, 
exp (— 2/c^/A^), has some energy dependence instead of being constant. It is to say, we replace 
in Eq. (2.7) G by G{E) = ci + C2k^, that reduces to the original form by setting C2 = 0. We 
find 

y^'-G = - + T4u-. + (-^ + A-.-T^]k' + ^^ + Oik^). (C.l) 

ci (27r)^/2 \^ c? CiA2 (27r)'^/2A/ 27r ^ ^ ^ ^ 

For model-given values of ci and C2 for an imposed cutoff A, one can shift the cutoff to A' and 
have the same T-matrix, up to 0{k^), by reabsorbing the cutoff shift in the new parameters c'l 
and c'2, given by: 

1 _ 1 , MA-AQ .p.^ 

c'l " ci (27r)3/2 ' ^^■^> 



c'2 _C2 , 2 A2 - A'2 2/i(A - A'^2 
' 7i A2A'2 ^ (27r)3/2AA' 



2 + ^ A2A/2 + /o_\.^rtA A,2 • (*^-3) 



If we insist in a constant potential, C2 = c'2 = 0, we can also reabsorb the cutoff effects in ci, 
but this would be correct just up to 0(/c^). 

Let us consider an approach in which the amplitude is written with a loop function regu- 



i*It satisfies [3], 



g2i<5 ^ 


k cot S + ik 
k cot 5 — ik 


Zo„{l,k^) + iTrik 


Zoo(l,fc2)-i7rifc 
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larized by means of a once-subtracted dispersion relation (DR) 
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^DR - ^DR ~ ^DR , (C.4) 

^DR = C , (C.5) 

Gdr = a -i-z- , (C.6) 

zir 

where C is the potential, analogous to the case of the Gaussian regulator approach, and a is a 
subtraction constant (a free parameter of the approach). Considering, as before, C = Ci +C2/c^, 
we can expand: 

^d"r - Gdr = - - « - ^A;2 + z/^ + Oik"") . (C.7) 



We can then reabsorb the effects up to 0{k^) of an arbitrary shift in the subtraction constant 
by means of: 

— = — + a'-a, (C.8) 

C'l Ci 

'2 

4 = %C2 . (C.9) 

Cl 

We see that the effects of the shift can be reabsorbed exactly for a constant potential, with 
the first of the previous equations. However, in the more general case of energy dependent 
potentials (as the case of chiral potentials, for example), this cannot be made exactly but just 
up to 0{k^). We see thus that there are several equivalent ways: one can fit, in the Gaussian 
regulator case, a constant for the potential and the cutoff, or fix the latter to a reasonable (but 
otherwise arbitrary) value and fit two constants. In a dispersion relation, a similar situation is 
found, where now the subtraction constant plays the equivalent role of the cutoff. In Subsec. 4.3 
we have followed the first approach. In Fig. 5, we show, for the case of / = J^'-' = O"*"^ DD, 
the contours curves of the x^ function in terms of the free parameters: the UV cutoff A and the 
inverse of the constant (Coa) that appears in the potential for this channel. We see already from 
this figure that A and 1/Coa are strongly correlated, and that the correlation is of the form 
given in Eq. (C.2). Indeed, the dashed line in the plot, that lies close to the axis of the error 
ellipse, is Eq. (C.2) using the central values of the cutoff and the potential given in the original 
work of Ref. [35], A = 1000 MeV and Coa = -1.024 fm^ We also infer that the correlation 
is stronger for higher values of the cutoff, since the quadratic terms in Eq. (C.l) become less 
important as A increases. 
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